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SUMMARY 


The response of dynamic systems to the product of two independent Gaussian 
random processes is developed by use of the Fokker-Planck and associated moment 
equations. The development is applied to the amplitude modulated process which 
is used to model atmospheric turbulence in aeronautical applications- The 
exact solution for the system response is compared with the solution obtained 
by the quasi-steady approximation which omits the dynamic properties of the 
random amplitude modulation. The quasi-steady approximation is valid as a 
limiting case of the exact solution for the dynamic response of linear systems 
to amplitude modulated processes. In the nonlimiting case the quasi-steady 
approximation can be invalid for dynamic systems with low damping. 


INTRODUCTION 

Random processes, which can be represented as the product of two inde- 
pendent random processes, or so-called product processes, include the ampli- 
tude modulated process which is used as a model for atmospheric turbulence in 
aeronautical applications (ref. 1). Atmospheric turbulence is considered to 
be a local Gaussian random process modulated by a relatively slow random ampli- 
tude process. Application of the amplitude modulated process requires methods 
for the analysis of the response of dynamic systems to the process. The cur- 
rently used method, termed the quasi-steady approximation in reference 2, con- 
siders the dynamic response of the system to the local Gaussian process only 
and assumes that the amplitude process affects the dynamic system only in a 
static manner. 

In the present study, the Fokker-Planck equation is used to develop a 
method of analysis that provides exact solution for the response of dynamic 
systems to the product of two independent Gaussian randcan processes. The method 
is used to obtain relations for the statistical moments of the response to the 
product process for some linear dynamic systems, including a linear time- 
invariant differential equation of general order. The exact and quasi-steady 
solutions for the response are conpared for both a first-order and a second- 
order system. The system response to the sum of the product process and an 
independent Gaussian process, proposed as an atmospheric turbulence model in 
references 3 to 6, is also examined. Some of the information in this report 
has been included in reference 7. 


SYMBOLS 

a linear filter constant of subscripted process 

b standard deviation of S process 

C( ) characteristic function 



e[ ] 

^jk 
f{ } 

f( ) 
hj 

h( ) 
i 

k,5, ,n 
M( ) 
M4 


m,M 


•N 


P( ) 
P(r| ro) 


q/Q 


r,R 


s,S 


t 

X 


y,Y 


z ,Z 


aj ( ) 


standard deviation of M process 
ensemble average 

response coefficient of linear state equation for system response 
Fourier transformation of indicated quantity 
response term of state equation 

excitation coefficient of state equation for system response 
= EjsrzjJ 

impulsive-response function of linear system 
unit of imaginaries, 
integer index (nonnegative) 
moment of indicated orders 

fourth-order flatness factor or kurtosis (fourth moment/square of 
second -moment) 

mean value randan process ' • 

sum of orders of moments' of system response 

probability density function 

conditional probability density function of r process conditional 
on value of ro 

uncoupled form of general local random process 
local randan process 
amplitude random process 
time 

general state vector random process 

response variable of second-order differential equation 
response variable of coupled linear system 

system response process; response variable of uncoupled linear system 
response term in general state equation 
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3j|^ excitation coefficient in general state equation 

Yn excitation coefficient in state equation for system response 

<S{t) Dirac delta function 

C Fourier transformation variable corresponding to z; system damping 

coefficient 

Dj Fourier transformation variable corresponding to yj 

X,v integer index (nonnegative) 

M linear filter constant for system response (negative of eigenvalue) 

Kj (t) Gaussian white noise process 

P Fourier transformation variable corresponding to r 

^ covariance matrix of coupled linear system 

o Fourier transformation variable corresponding to s 

T time-difference variable 

w system frequency coefficient 

Subscripts: 

in input 

j,k,8-,m integer index 

o at an earlier value of time 

out output 

The summation convention is used; repeated indices imply a summation 
unless noted otherwise. 

A dot over a symbol indicates a derivative with respect to time. 

RESPONSE OF DYNAMIC SYSTEMS TO PRODUCT PROCESSES 

Product randcm processes, particularly the special case of the amplitude 
modulated process, are described in this section. The quasi-steady approxima- 
tion for the analysis of the response of dynamic systems to the amplitude modu- 
lated random process is examined. The quasi-steady approximation, which 
ignores the dynamic properties of the amplitude process, is described in quali- 
tative terms. The Fo)<ker-Planck equation is developed for the response of a 
first-order nonlinear dynamic system to the product process. The response of 
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the system to the sum of a product process and an independent mean value pro- 
cess is also considered. 


Product Processes and the Quasi-Steady Approximation 

The product randan process is formed by the product of two canponent 
processes 


r(t) s(t) (A) 

The component processes are specified to be stationary, independent, and 
Gaussian with zero mean values. The general product process (expression (A)) 
includes the special case of the amplitude modulated process where the local 
R canponent is a rapidly varying process which is modulated by the slowly 
varying amplitude component S.^ The product process is developed in general 
terms without restriction on the relative time variations of the two canponent 
processes. Since the probabilistic structures of the two components are com- 
pletely defined (given their autocovariance functions), the corresponding struc- 
ture of the product process is determined by the defining relation (expres- 
sion (A)) and by the appropriate relations of probability theory (ref. 2). 

The amplitude modulated process was originally developed by consideration 
of measured atmospheric turbulence data, segments of which appear to be sta- 
tionary and Gaussian over small time intervals. Over longer time intervals, 
the amplitude of the local Gaussian segments appears to vary in a randan manner. 
This interpretation of measured turbulence data led to the development of the 
amplitude modulated process, which has had extensive application in aeronautical 
problems associated with atmospheric turbulence (refs. 8 to 10). 

The amplitude modulated randan process has been applied in the analysis 
of the response of dynamic systems. The associated analytical technique, termed 
the quasi-steady approximation, uses the idea of the relatively slow variation 
of the amplitude component to develop a simple method of analysis. The develop- 
ment of the quasi-steady approximation is shown by consideration of the response 
or output of a linear system 

zout(t) = h(T) Zin(t-T) dT (1) 

The system response is the convolution of the impulsive response function of 
the linear system and the excitation or input function, which is a product 
process (expression (A)) 


^The terms product process, amplitude modulated process, and quasi-steady 
process have the following meanings; a product process includes all randcsn 
processes of the form of expression (A) ; an amplitude modulated process is the 
special case with the S canponent slowly varying with respect to the R 
component process; a quasi-steady process is the limiting case where the 
dynamic properties of the S canponent are omitted. 
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^OUt (t) 


h(T) rin(t-T) Sin(t-T) dx 


(2) 


= 1 


Two assumptions are made. First, it is assumed that the impulsive response 
function is approximately equal to zero after a known time period. Conse- 
quently, the integration in equation (2) can be restricted to this period with 
a negligible change in the computed system response. Second, it is assumed 
that the amplitude component process is approximately constant during this time 
period. With these two assumptions the amplitude process can be removed from 
the convolution integral to give 


^out “ ®in (t) \ h(T) rin(t-T) dT (3) 

The remaining integral is the response of the linear system to the local 
R component process. Since the input R process is Gaussian, the output R 
process is also Gaussian. The approximate system response satisfies the pro- 
duct relation and consequently is an amplitude modulated process 




(4) 


The dynamic aspects of the response result from the rapidly varying R compo- 
nent process only. The amplitude component influences the system only in a 
static manner. 

The preceding development shows the conditions that are required for the 
quasi-steady approximation. The amplitude component must be slowly varying 
relative to both the local R component and the dynamic properties of the 
system. The accuracy of the approximation is improved if the system is well 
damped; the accuracy is questionable if the system is lightly damped. The 
quasi-steady relation is exact for static linear systems; this follows from a 
linear transformation of the defining relation (expression (A)). The quasi- 
steady method of analysis is generally an approximation, however, since devel- 
opment of the exact response requires the consideration of the dynamic response 
of the system to both component processes. If the exact response is known, 
then the quasi-steady response is the limiting case as the dynamic properties 
of the amplitude process become negligible. 

An alternate random process has been proposed as a model of atmospheric 
turbulence in references 3 to 6. The process is formed by the sum of a product 
process and an independent, random variation of the mean value 


r(t) s(t) + m(t) 


(B) 


The resulting product-plus-mean process is formed from three component processes 
which are specified to be stationary, independent, and Gaussian with zero mean 


5 



values. Since the probabilistic structures of the three component processes 
are completely defined, the structure of the product-plus-mean process is deter- 
mined by the defining relation (expression (B) ) and the appropriate relations 
of probability theory. 

The response of dynamic systems to the product-plus-mean process can be 
determined by the quasi-steady approximation. The addition of the mean value 
component presents no fundamental difficulty for the analysis of linear dynamic 
systems, since the response to the independent product and mean value processes 
can be considered separately. The analysis of the response to the Gaussian 
mean value process requires only the determination of the covariance functions 
of the response. Using the quasi-steady approximation for the response to the 
product process, the system response is also a product-plus-mean process 
(expression (B) ) . For nonlinear dynamic systems the analysis is inherently 
more difficult, since the response of the system to the product and the mean 
value processes cannot be considered separately. 


Pokker-Planck Equation for the Response of Dynamic Systems 

The general Fokker-Planck equation is applied to the response of dynamic 
systems forced by the product process. The response of a first-order nonlinear 
dynamic system is considered in order to illustrate the general procedure. The 
development is related to that in reference 11. 

In the present example the two components of the product process are gen- 
erated by first-order filtering of independent Gaussian white noise processes. 
The corresponding differential equations are 


s. = -agS + ^^2^ b K-\ (t) (5a) 

r = -aj-r + ^^2a7 ? 2 (t) (5b) 


The two filter constants as and a^ determine the dynamic properties of the 
two component processes. The two independent Gaussian white noise processes 
have zero mean values and have autocovariance functions which are Dirac delta 
functions 


o 

II 

1 


-1 

[5k(tl) C£(t2)J =< 

(k = H) 
(k fi 1) 




( 6 ) 
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The corresponding power spectral density functions of the white noise processes 
have a constant value of 1/2 tt. The variances of the two conponent processes 
are 

e[s2J = 

E[r2] = 1 J 

A first-order nonlinear dynamic system with excitation by the product 
randan process (expression (A) ) is considered. 


z = -f (z) + rs 


( 8 ) 


A composite dynamic system is formed which consists of the differential 
equations for the dynamic system (eq. (8)) and for the conponent processes 
(eqs. (5a) and (5b)). These equations are combined to form the nonlinear state 
equation for the composite system, which has the general form 


Vj = ttj (V) + Bjg, CjJ,(t) 


(9) 


The state vector v is the array of the variables of the composite system, 
which in the present case includes the two components of the excitation process 
and the system response, that is, the three variables s, r, and z. The 
excitation term of the composite dynamic system consists of the two Gaussian 
white noise processes of equations (5a) and (5b). The excitation term of the 
original dynamic system (eq. (8) ) is the product of two response variables of 
the composite system and, consequently, is included in the response term, that 
is, the a functions of the state equation (eq. (9)). 

The probability density function of the composite system satisfies the 
associated Fokker-Planck equation (refs. 12 to 15). For a dynamic system satis- 
fying the general state equation (eq. (9)), the Fokker-Planck equation is 


3t 




a2 

9vj 9v|^ 




( 10 ) 


The Fokker-Planck equation determines the transition probability density func- 
tion of the state variables. 





) 
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which is the probability density function of the state vector conditional on 
the given value of the state vector at an earlier value of time. The transition 
probability density function is generally nonstationary due to the development 
of the state vector frcm its initial value. ^ if the conditional dependence on 
the initial value is omitted, the resulting first-order probability density 
function is stationary; the time derivative term in the Fokker-Planck equation 
is equal to zero. 

The general Fokker-Planck equation (eq. (10)) is applied to the analysis 
of the response of the composite dynamic system to the product process. The 
associated probability density function is the conditional joint density of the 
two component processes and the system response 

p(s,r,z| So,ro,Zo) 

The Fokker-Planck equation for the composite dynamic system of equations (5a) , 
(5b), and (8) is obtained by application of equations (9) and (10). The result- 
ing equation is 


3t 


3 , 3^P 


3 3 2p 

— (sp) + b^ 

+ 

— (rp) + 

_3s 3 s2_ 


_3r 3r2_ 



9p 
rs — 
3z 


( 11 ) 


The first two terms of the right-hand side correspond to the differential equa- 
tions of the S and R component processes. The next term corresponds to the 
nonlinear response term of the original dynamic system (eq. (8)). The final 
term corresponds to the excitation of the original dynamic system by the product 
process. This term couples the system response and the two components of the 
excitation process in the Fokker-Planck equation. 

A similar procedure is used for the analysis of the response of dynamic 
systems to the product-plus-mean process. An additional equation defining the 
mean value process is introduced. In the present example, the mean value pro- 
cess is generated by a first-order filtering of a Gaussian white noise process 
in the same manner as the two other component processes 


m = -a„,m + ^2a„, c ^3(t) 


( 12 ) 


^The development of a random process frcm its initial value introduces the 
concept of transition properties. These properties are particularly important 
for the amplitude modulated process due to the significant difference in the 
dynamic properties of the two component processes. The transition properties 
explicitly show that short samples of the process appear to be Gaussian with 
an approximately constant amplitude factor and that long samples have the form 
of the product process (ref. 2) . 
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The differential equation for the system response is modified to account for 
excitation by the combination of the product and the mean value processes 


z 


-f (z) + rs + m 


( 13 ) 


The Fokker-Planck equation is developed from the general form (eq. (9)) 
of the state equation for the composite dynamic system. The associated proba- 
bility density function is the conditional joint density of the three component 
processes and the system response 


p(sfr,m,z SQ,rQ,mQ,ZQ^ 


The Fokker-Planck equation for the composite dynamic system of equations (5a), 
5(b) , (12) , and (13) is 


3t 


as 


3 

— (sp) 
3s 



+ 


ar 


3 


3r2 


+ 


^m 


3 

— (mp) 
3m 




3p 3p 

rs — - m — 

3z 3z 


(14) 


The resulting Fokker-Planck equation is similar to equation (11) except for the 
presence of two additional terms; one term corresponds to the differential 
equation for the mean value process and one term corresponds to the additional 
excitation of the system by that process. 


RESPONSE OF LINEAR DYNAMIC SYSTEMS 

Several examples of the response of linear dynamic systems to the product 
random process are developed in the present section. Exact relations for the 
statistical moments of the system response are obtained by using the Fokker- 
Planck equation for the associated composite dynamic system. The quasi-steady 
approximation is developed as a limiting case of the exact solution for the 
system response. The analysis shows the conditions required for the approxima- 
tion in the dynamic systems considered. 

The method is first applied to the analysis of the response of a first- 
order linear dynamic system to the product process. The method is then used 
to develop the solution for a linear dynamic system of general order, with 
application to a second-order linear differential equation. One example of the 
response of a linear dynamic system to the product-plus-mean process is con- 
sidered. Finally, several possible extensions of the development are outlined. 
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Response of First-Order Linear System to Product Process 


The response of a first-order linear system to the product process, with 
both components generated by first-order filtering of white noise processes, 
is considered. The composite dynamic system consists of three first-order dif- 
ferential equations; two are for the component processes and one is for the 
system response as follows: 


s = -agS + \j2ag b ?-| (t) 

(15a) 

r = -a^r + \j2aj- ^2 (t) 

(15b) 

« 

z = -pz + rs 

(15c) 


The three filter constants (ag, a^^, and y) determine the dynamic properties 
of the two components of the product process and those of the system response. 
The filter constants are positive and are the reciprocals of the time constants 
of the associated first-order filters. With this system, the dynamic properties 
of the amplitude process can be varied relative to those of both the R compo- 
nent process and the system response by changing the values of the filter 
constants. 


The frequency response functions of the three state variables of the compo- 
site dynamic system are shown in schematic form in figure 1 . Since the two 
excitation components both satisfy first-order linear differential equations 
with constant coefficients, their frequency response functions have an approxi- 
mately constant value for frequency values below the corresponding filter con- 
stants and decrease as a -1 power of the frequency above that value. The fre- 
quency response function of the system response has a similar form.^ The 
frequency response function of the first derivative of the system response 
varies as a +1 power of the frequency below the value of the filter constant 
and is constant above that value of the frequency. The frequency response func- 
tions in figure 1 are shown for the modulation condition with the amplitude 
process S slowly varying relative to the local R process, that is, 
as « aj-. 

Using equations (9) and (10), the Fokker-Planck equation for the composite 
system of equations (15a) to (15c) is 


9p 




+ 


ar 


3 

— (rp) 
3r 


+ 


32p 


3 

+ M — (zp) 
3z 


3p 
rs — 
3z 


(16) 


^Figure 1 shows the frequency response function for the first-order system 
of equation (15c), that is, for the system response relative to the excitation 
product process. 


10 



The probability density function is the conditional joint density of the three 
state variables: s, r, and 2 . Equation (16) is a special case of equa- 

tion (11). The first two terms on the right-hand side correspond to the dif- 
ferential equations for the tv/o ccmponent processes (eqs, (15a) and (15b)). 

The last two terms correspond to the differential equation of the system 
response (eq. (15c)), the last term corresponding to the excitation of the 
system by the product process. 

The Fokker-Planck equation is transformed to give the corresponding equa- 
tion for the joint characteristic function of the composite dynamic system. 

The characteristic function is the Fourier transformation of the probability 
density function over the three state variables 


C P, ? 1 Sq, r Q, Zq^ — F^p^s,r,z 


®0' O' 





(17) 


Taking the Fourier transformation of the Fokker-Planck equation, the resulting 
characteristic function equation is 


9C 

3 g 

r p 5 1 

a — + b^a-^C 

+ a^ 

3C 

P — + P'^c 

9t ~ 


9a 


3p 


9C a 2c 

+ r-T" 

9c 3a 9p 


(18) 


(The development of this equation requires relations which are derived in 
appendix A.) The terms of the characteristic function equation correspond to 
the terms of the Fokker-Planck equation. The final term corresponds to the 
excitation of the dynamic system by the product process. 

Relations for the statistical moments of the system response can be devel- 
oped frcm either the Fokker-Planck or the characteristic function equation. 

The following notation is used for the joint moments of the three state 
variables: 


E[s*^r^za]] = M(k,S-,n) 


(19) 


The moment equation is obtained frcm the characteristic function equation by 
taking the derivatives of general order with respect to the Fourier transforma- 
tion variables, which are then set equal to zero to develop the moments from 
the characteristic function. Using this procedure and some relations developed 
in appendix A, the equation for the joint moments is 
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( 20 ) 


3 

— M(k,S, ,n) 
3t 


= k(k - 1) 


+ (kag + Zaj- + nM ) M(k,(!-,n) 
agb^M (k-2 ,2. ,n) + 2. (2< - 1) a^-M (k,2,-2,n) 


+ nM (k+1 ,2/+l ,n-l ) 


The individual terms of the moment equation correspond to the individual terms 
of the state equation for the composite dynamic system. The first term on the 
left-hand side corresponds to the transient response of the ccmposite dynamic 
system. The other terms on the left-hand side correspond to the three linear 
response terms. The first two terms on the right-hand side correspond to the 
two white noise excitation terms. These terms pass the moments of the S and 
R processes from lower to higher order without changing the order of the moment 
of the system response. The final term corresponds to the product excitation 
term in the differential equation for the system response. This term introduces 
coupling between the moments of different orders in the system response. 


The moments of the stationary response of the systoti are obtained by 
deleting the time-derivative term in the moment equation. The subsequent dis- 
cussion is restricted to the stationary response. Examination of the moment 
equation shows that the odd-order stationary moments of the system response are 
zero. The second moment (which equals the variance in this case) of the system 
response is 


E z2 = (21) 

M (ag + aj. + M) 

The variance of the derivative of the system response is obtained from equa- 
tions (15c) and (21) 


E 



(ag + ai-)b2 
(ag + aj- + M) 


( 22 ) 


The variances of the system response and its derivative depend upon the filter 
constant ag of the amplitude process. If the value of this filter constant 
is much smaller than that of the R process, then the variances become those 
obtained by the quasi-steady approximation. The quasi-steady approximation 
gives an accurate prediction of the variances of the system response and its 
first derivative in the limiting case (ag/a^ = 0) irrespective of the value of 
the filter constant )j of the dynamic system. 

The higher moments of the system response are computed from the moment 
equation (eq. (20) ) by a procedure which is described subsequently. In fig- 
ure 2, the resulting fourth moments of the system response and its first deriva- 
tive are plotted in normalized form, that is, normalized to their quasi-steady 
values. The moments are plotted as functions of the filter constant of the 
amplitude conponent and for several values of the filter constant of the linear 
system, both constants given as ratios to the filter constant for the local 
R component process. The quasi-steady approximation corresponds to the zero 
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value of the abscissa. The dynamic properties of the amplitude component 
reduce the value of the fourth moment of the system response. The fourth 
moment of the derivative can be either decreased or increased by the dynamic 
properties of the amplitude ccmponent. For both the system response and its 
first derivative the value of the fourth moment is decreased as the value of 
the filter constant of the dynamic system is decreased. 

The moments of the system response are also presented in terms of the 
associated flatness factors, which indicate the relative occurrence of the 
extreme values of the random process. The values of the fourth-order flatness 
factor (or kurtosis) of the system response and its first derivative are plotted 
in figures 3 and 4. The flatness factors are plotted in the same manner as fig- 
ure 2, that is, as functions of the filter constant of the amplitude process 
and for several values of the filter constant of the dynamic system, both con- 
stants given as ratios to that of the local R process. In the quasi-steady 
limit, the values of the two flatness factors both approach the value of 9 for 
the product of two independent Gaussian processes, which is the quasi-steady 
approximation (eq. (4) ) . The values of the flatness factors are reduced by the 
dynamics of the amplitude process. 

The results of figures 2 to 4 indicate whether use of the quasi-steady 
approximation will either overestimate or underestimate the occurrence of the 
extreme values of the system response and its first derivative. Since the 
values of the fourth-order flatness factors are maximum in the quasi-steady 
limit, use of the quasi-steady approximation in combination with the exact 
values of the second moments (eqs. (21) and (22)) will tend to overestimate the 
occurrence of the extreme values. If the quasi-steady values of the second 
moments are used, then the values of the fourth moments (fig. 2) and not those 
of the fourth-order flatness factors must be examined. For the system response, 
use of the quasi-steady approximation in this form will tend to overestimate 
occurrence of the extreme values. For the derivative of the system response, 
use of the quasi-steady approximation can either overestimate or underestimate 
the occurrence of the extreme values, depending upon the value of the ratio of 
the filter constants of the first-order system and the R process. 

The exact solution for the response of the first-order system to the pro- 
duct process is interpreted in the following manner: For the system response 

(fig. 3) , the value of the fourth-order flatness factor is increased and 
approaches the quasi-steady value of 9 as the value of the system filter con- 
stant M is increased. In the limit of large values of the filter constant, 
the system response becomes static with respect to both components of the 
product process, since the frequency response function is approximately con- 
stant for frequency values below n (fig. 1). Consequently, the system 
response is approximately a product process resulting from a static transforma- 
tion of the excitation product process. For the derivative of the system 
response (fig. 4) , the trend is reversed. The value of the fourth-order flat- 
ness factor increases and approaches the quasi-steady value of 9 as the value of 
y is decreased. In this limit, the derivative becomes a static transformation 
of the product process and is thus a product process. This limiting property 
of the derivative also follows directly from the differential equation of the 
original dynamic system, that is, by setting the value of y to zero in 
equation (15c). 
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The stationary forms of the Fokker-Planck and characteristic function equa- 
tions (eqs. (16) and (18)) can be solved in two special cases. In the first 
case, the value of the parameter C in the characteristic function equation is 
set to zero. The resulting equation is that for the characteristic function of 
the two independent Gaussian processes S and R. In the second case, the 
value of the filter constant ag of the amplitude process is set to zero in 
the Fokker-Planck equation to obtain the special form for the quasi-steady 
approximation, which is 


a 


a 


r 


— (rp) 
3r 


+ 


a2p 


a 

+ y — (zp) 

az 


3p 

rs — = 
3z 


0 


(23) 


The probability density function is replaced by the conditional function, which 
is the joint probability density function of the system response and the R 
component, conditional on the value of the amplitude parameter s. This condi- 
tional joint probability density has the Gaussian functional form with the 
appropriate covariance matrix for the system response and the R component. 


Response of General Linear System to Product Process 

The relations for the moments of the response of a linear time-invariant 
dynamic system of general order are now developed. Relations are developed 
first for the response of a set of uncoupled first-order systems, the develop- 
ment being a direct extension of that for the single first-order system con- 
sidered previously. The relations for the uncoupled system can be applied to 
a general linear system by introducing the eigenvalues and eigenvectors of 
the coefficient matrix of the corresponding state equation to uncouple the 
general linear system. The response moments for the general system can thus 
be ccanputed from those for the uncoupled system. Relations for the covariance 
matrices of the response of the general linear system are also developed. 

These relations can be expressed without the use of the eigenvalues and eigen- 
vectors of the state coefficient matrix. 


Moment equation .- For a set of uncoupled first-order linear systems, the 


component 

processes and 

the equation for the uncoupled system. 


= -ags + 

^2as b Cl (t) 


(24a) 

= -aj-r + 

\l2aj. C2(t) 


(24b) 

j = -UjZj 

+ rsYj 

(No summation on j) 

(24c) 
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This composite system is a direct extension of equations (15a) to (15c). The 
associated filter constants, which are the negatives of the eigenvalues, can 
be complex with the restriction that the real parts be positive. The moment 
equation for the composite system is obtained by developing the associated 
Fokksr "Planck and characteristic function equations. The development is similar 
to that of equations (16) to (18). The equation for the characteristic function 
of the composite system of equations (24a) to (24c) is 


3C 

^s 

ac 

a — + b^a^c 

+ aj. 

' 2 ' 
p — + p-^c 

at ' 


aa 


L 9p J 


3c a^c 

aci ao ap 


(25) 


The Cj are the Fourier transformation variables corresponding to the Zj vari- 
ables. The corresponding moment equation is 


— + kag + Siaj. + ^ VjVij 

j=l 


M(k;fi,;V-] 


•»Vn) 


= k(k - 1) agb2 M(k-2;5,;Ui 


V„) + jl(Jl - 1) ai-M(k;S,-2;vi , . . . , Vn) 



VjYj M(k+1 ;2,+l ;Vi 


(Vj-1) 


• • * 


>^n) 


(26) 


where 


M(k;H;V] ,V2, . . . /V^) 


e[s^ r^ 



The individual terms of the moment equation correspond to the individual terms 
of the state equation for the composite system. The terms on the left-hand 
side, aside from the time-derivative term, correspond to the linear response 
terms in the set of first-order differential equations. The first two terms 
on the right-hand side correspond to the two white noise processes. The final 
set of terms on the right-hand side corresponds to the excitation of the 
uncoupled dynamic system by the product process. 

In the stationary case, the time-derivative term of the moment equation 
is zero. The stationary moments of the system response can be computed from 
equation (26) by a procedure which is based upon the sum of the orders of the 
moments of the original dynamic system 


n 
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The moments are developed by computing all combinations of moments for a given 
value of the index N, and then increasing the value of that index by one. The 
procedure starts with the index N having the value zero; the corresponding 
moments are those of the independent and Gaussian R and S component pro- 
cesses, As the value of the index N is increased, the moments are generated 
by the first two terms on the right-hand side of equation (26) (which introduce 
the moments of lower order in the R and S components) and by the final set 
of terms on the right-hand side (which passes the moments from one value lower 
in index N (and therefore known) and from one value higher in the indices of 
the R and S components) . This procedure is easily implemented on a digital 
computer. The only special consideration is the core storage required if the 
size (index n) of the dynamic system is large. 

In the nonstationary case, the moment equation is used to compute the time- 
dependent transition moments which result from the development of the system 
response from a given initial condition. Given the moments at one time value, 
the moment equation is used to conpute the first derivatives (with respect to 
time) of the moments. The derivatives are used to determine the moments at a 
subsequent time value. The process is then repeated. 

The moment equation of the set of uncoupled first-order systems (eq. (24c)) 
can be used to conpute the moments of a coupled linear system 


yj = FjkVk + rsgj 


(28) 


The transformation between the two systems is accomplished by using the eigen- 
values and eigenvectors of the state coefficient matrix F of the coupled sys- 
tem, The moment equation is thus the essential element in the development of 
the moments of the response of a general time-invariant linear system to excita- 
tion by the product of two Gaussian processes. 

Covariance matrix equation .- Specific relations can be developed for the 
covariance matrices of the response of a general linear system. Since the 
equation has a reasonably simple form, the covariance matrices can be computed 
directly without the introduction of the eigenvalues and eigenvectors of the 
system. The development of the covariance matrix equation is a standard method 
for the analysis of the response of linear dynamic systems to Gaussian white 
noise processes (refs. 14, 16, and 17). 

The state equation of the composite dynamic system is formed from equa- 
tions (24a), (24b), and (28). The corresponding Fokker-Planck , characteristic 
function, and moment equations are developed in the same manner as for the 
first-order systems considered previously. The joint probability density func- 
tion of the system response and the two component processes is determined from 
the Fokker-Planck equation for the conposite dynamic system 


— 

3t 


3 3V 

(sp) + b2 

3s 3 s2 


3 32p 

+ ar — (rp) + 

3r 3r2 




(29) 
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The joint characteristic function is the Fourier transformation over all inde- 
pendent variables of the joint probability density function 


c(u,P,r,j |so,ro,yjo) 


= i^(p(s,r,yj |so,ro,yjo)} 


The corresponding characteristic function equation is the Fourier transformation 
of the Fokker-Planck equation 



The equation for the covariance matrix of the system response variables 
is obtained by taking the partial derivatives of the characteristic function 
equation with respect to two system variables and then developing the moments 
from the characteristic function 


3 

- — ^j£ + + Ffi,kZkj = -gjhj, - giihj (32) 

where 

^kS- = E[yky£] 
hj = E[sry3] 

This equation is similar to that for the response of linear dynamic systems to 
white noise processes and differs only in the form of the term on the right- 
hand side, which is related to the excitation process. The solution of the 
equation for the covariance matrix of the response requires the array of triple 
moments of the two conponent processes and the array of state variables of the 
original dynamic system. The equation for these moments is obtained from the 
characteristic function equation by taking the partial derivatives with respect 
to a, p, and one of the n variables and then developing the moments from 
the characteristic function 

a „ 

hj + (ag + aj-)hj - Fj^h^ — b gj (33) 

This equation can be solved directly for the array of triple moments, which 
is then used to determine the covariance matrix of the system response by 
equation (32) . 
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The covariance matrix of the first derivative of the system response is 
obtained fron the state equation of the system response (eq. (28)) 


E[yjk] = Fj£^ilmFkm + FjiLhjigk + Fkj;,h£gj + b2gjg,^ (34) 

This covariance matrix can be determined once that of the system response is 
known. 


The equations for the covariance matrix (eqs. (32) and (33)) show the 
effect of the quasi-steady approximation. The filter constants of the amplitude 
and local R processes appear only in their sum in equation (33). Thus, if 
the amplitude process is much slower than the R process, that is, ag « a^^, 
the error in using the quasi-steady approximation (ag/a^ = 0) is small for the 
covariance matrix. Also, the quasi-steady approximation does not depend on the 
dynamic properties of the linear dynamic system in the limiting case, although 
those properties generally affect the accuracy of the quasi-steady approximation 
in the nonlimiting case- 


Stationary Response of Second-Order Linear System 

The preceding development is applied to the analysis of the response of 
a second-order linear system, to the product process. The differential equation 
for the dynamic system is 


X + 250JX + 


0)2x = 


rs 


(35) 


The state equation for the composite dynamic system is developed by writing 
equation (35) as a state equation and combining with the first-order equations 
for the S and R conponents (eqs. (24a) and (24b) ) . The covariance matrix 
equation is then written for the composite dynamic system by using equa- 
tions (32) and (33). The covariance matrix equation is solved for the vari- 
ances of the response and its first derivative; this gives 



b2 




Bto2L 


(ag + aj.) 
+ 


E 



(ag + ar)b2 

2?(jj3 


(36a) 


(36b) 


where 3 = (ag + aj-)2 + 2Cw(ag + a^.) + 0)2. The quasi-steady approximation is 
obtained from the exact solution by setting the ratio ag/a^. to zero. The 
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quasi-steady form of the two variances may be either larger or smaller than 
the exact values, depending on the values of the system damping and frequency 
coefficients. 

The higher moments of the system response are developed from the general 
relation (eq. (26)) by using the eigenvalues and eigenvectors of the state 
equation of the second-order linear system. In figures 5 and 6, the resulting 
fourth moments of the system response and its first derivative are plotted in 
normalized form, that is, normalized to their quasi-steady values. The moments 
are plotted as functions of the filter constant of the amplitude process and 
for several values of the system frequency coefficient to, both quantities 
plotted as ratios to the filter constant of the local R process. The moments 
are plotted for two values of the system damping coefficient C. The values 
of the fourth moment are usually decreased by the dynamic properties of the 
amplitude process. The relative decrease is larger for smaller values of the 
frequency coefficient to. In the case of large values of the ratio (o/aj- , the 
dynamic properties of the amplitude process can increase the moments of the 
response and particularly the first derivative of the system response. The 
effects of the dynamic properties of the amplitude process are usually more 
significant for lower values of the system damping coefficient C- 

The fourth-order flatness factors of the system response and its first 
derivative are shown in figures 7 to 12. These are plotted in the same manner 
as the fourth moments in figures 5 and 6. Figures 7 and 8 show the flatness 
factors of the system response and its first derivative for the system damping 
coefficient C of 0.707. In both cases the fourth-order flatness factor has 
the value of 9 in the quasi-steady limit (ag/aj- = 0) . The values of the flat- 
ness factors are reduced by the dynamic properties of the amplitude process. 

For the system response (fig. 7) , the value of the flatness factor is decreased 
as the value of the system frequency coefficient oo is decreased. As the value 

of the system frequency is increased the flatness factor approaches the quasi- 
steady value of 9. The qualitative features of the response of the second-order 

system are thus similar to those of the response of the first-order system 

(fig. 3). The flatness factors of the first derivative of the system response 
(fig. 8) show the same features as the response itself except that the flatness 
factor is less sensitive to the value of the system frequency coefficient w. 
Also, the flatness factor of the first derivative does not approach the quasi- 
steady value of 9 in the limit of large values of the system frequency coeffi- 
cient. Figures 9 and 10 show the fourth-order flatness factors of the system 
response and its first derivative for the system damping coefficient C of 0.10. 
The values of the flatness factor tend to be considerably less than those for 
the higher value of the damping coefficient in figures 7 and 8. Otherwise the 
flatness factors show the same general features, including the value of 9 in 
the quasi-steady limit. 

The results of figures 5 to 10 indicate whether use of the quasi-steady 
approximation will either overestimate or underestimate the occurrence of the 
extreme values of the system response and its first derivative. Since the 
values of the fourth-order flatness factors are maximum in the quasi-steady 
limit, use of the quasi-steady approximation in combination with the exact 
values of the second moments (eqs. (36a) and (36b)) will tend to overestimate 
the occurrence of the extreme values of both the system response and its first 
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derivative. If the quasi-steady values of the second moments are used, then 
the quasi-steady approximation will tend to overestimate the occurrence of 
extreme values of the system response (fig. 5) . (In the limit of large values 
of the ratio w/aj- the quasi-steady relation becomes exact, since the system 
response becomes a static transformation of the excitation product process.) 

The same form of the quasi-steady approximation can either overestimate or 
underestimate the occurrence of extreme values of the first derivative, depend- 
ing upon the value of the ratio w/aj. (fig. 6). 

The effect of the damping coefficient is shown explicitly in figures 11 
and 12 where the fourth-order flatness factors are plotted for several values 
of the damping coefficient C and for a single value of the system frequency 
coefficient. For both the system response and its first derivative, the flat- 
ness factor is strongly dependent upon the value of the damping coefficient. 

As the value of the damping coefficient is increased, the flatness factors are 
increased except for values of the damping coefficient near 1. For low values 
of the damping coefficient, the flatness factors of both the response and its 
first derivative are close to the Gaussian value of 3. The improvement of the 
quasi-steady approximation with increased system damping was indicated by the 
qualitative development of the quasi-steady approximation (eq. (4)). In all 
cases the appropriate value of 9 is obtained in the quasi-steady limit. For 
small values of the system damping coefficient, the dynamic properties of the 
amplitude process invalidate the quasi-steady approximation, that is, for small 
(but nonzero) values of the ratio of the filter constants ag/aj. the flatness 
factors approach the value of 3. Thus for lightly damped second-order linear 
systems, the dynamic properties of the amplitude component process largely 
destroy the effect of the amplitude modulation. 

A general rule of random process theory is that linear filtering of a non- 
Gaussian stationary process results in a tendency toward a Gaussian stationary 
process (refs. 18 and 19). This rule is based upon the expression of the linear 
system response, using the impulsive response function and the convolution 
integral, as a weighted sum of the non-Gaussian excitation variables over a 
range of time values. Since the summation introduces an element of statistical 
disorder, there is a tendency of the sum toward a Gaussian distribution as indi- 
cated by the central limit theorem. Further, this tendency is stronger if the 
system response has a longer duration, that is, the system frequency is lower 
or particularly if the system damping is lower. The previous results for the 
fourth-order flatness factors for the response of the second-order system to 
the product process are consistent with these rules. 


Response of First-Order System to Product-Plus-Mean Process 

An example of the analysis of the response of a linear dynamic system to 
the product-plus-mean process (expression (B) ) is briefly considered. Since 
the dynamic system is linear, the response to the independent product and mean 
value processes can be considered separately. Since the mean value process is 
Gaussian, the corresponding linear system response is also Gaussian. This 
response is added to that for the product random process. 
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The linear system of equation (15c), modified by the addition of the mean 
value component into the excitation process, is considered for an example. The 
mean value process is formed by a first-order filtering of an independent white 
nol.so pirocsss in ths iriwnnojr 3S ^1*0 Xqo 3!L R cGnponsnts 


m 


-amm + c ^3<t) 


(37) 


The dynamic system is a first-order linear system excited by the product-plus- 
mean process 


-JJz + rs + m 


(38) 


The analysis of the system response can be developed by using the Fokker- 
' 'Planck, characteristic function, and moment equations for the composite dynamic 
system. However, it is easier to consider the response to the independent pro- 
duct and mean value processes separately 


z = zi + Z2 


(39) 


The first term on the right-hand side is the response to the product process; 
the moments are given by equations (21 ) and (22) . The second term is the 
response to the mean value process. The second moment (which equals the vari- 
ance in this case) of the total response is 


E 



b2 

y (ag + aj- + M) 



y (ajn + y) 


(40) 


The covariance of the derivative of the response is found by using the differen- 
tial equation of the original system (eq. (38)) and equation (22) 


E 



(ag + aj.)b2 



+ y) 


3mC^ 

(am + y) 


(41) 


These two relations show the relative influence of the product and the mean 
value processes,/ which are identified by their standard deviations b and c, 
respectively. All other moments of the system response can be developed in a 
similar manner./ The main problem is the analysis of the system response to the 
product process; the analysis of the additional response to the mean value pro- 
cess introduces no fundamental difficulty for linear systems. 
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Possible Extensions of the Analysis 


The present development of the exact analysis of the response of dynamic 
systems to product random processes can be extended in several ways. One exten- 
sion is the development of the analysis for more general forms of the two compo- 
nent processes. The primary requirement is that the components be generated 
by linear time-invariant filtering of Gaussian white noise processes. The dif- 
ferential equations for the components become more complicated, but the formula- 
tion procedure for the associated Fokker-Planck, characteristic function, and 
moment equations is not changed. The relations for the response of a linear 
system to a product process having one component formed by a general-order 
filtering of Gaussian white noise are developed in appendix B. 

The present development can be extended to the analysis of dynamic system 
response to more general forms of the product process. One possible form, for 
example, is the sum of several independent Gaussian and product processes 


“I = 


mj (t) 



Sj (t) 


All component processes are specified to be stationary, Gaussian, and indepen- 
dent. The general form of the preceding expression allows considerable flexi- 
bility in modeling both the dynamic and probabilistic structures of a random 
process. The response of a linear dynamic system to the general random process 
is the sum of the terms for the response to the individual terms of the excita- 
tion process, where the individual terms of both the excitation and response 
processes are independent. The moments of the system response to the individual 
terms and thus to the total excitation process can be determined by the present 
analytical method. 

The present development can also be extended to the analysis of the non- 
stationary response of dynamic systems to product processes. The moments are 
functions of time due to the nonstationarity of the system response. One 
example of the nonstationary case is the transient response of linear time- 
invariant systems due to initial conditions (refs. 14 and 16). A second 
example is the analysis of the response of linear dynamic systems with time 
dependent coefficients (ref. 17). 


CONCLUDING REMARKS 

Exact and approximate solutions for the response of dynamic systems to the 
product of two independent and stationary Gaussian processes are examined. The 
Fokker-Planck equation is used to develop the exact solution for the system 
response. For linear time-invariant systems, the associated moment equation 
can be formulated and solved for all statistical moments of the system response. 
The development thus presents a method for the analysis of the response of 
dynamic systems to a class of strongly non-Gaussian random processes. 
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The validity of the quasi-steady approximation for. the analysis of the 
response of dynamic systems to amplitude modulated processes is examined. The 
quasi-steady solution is developed from the exact solution by omitting the 
dynamic properties of the modulating amplitude randcm process. The quasi- 
steady approximation requires that the amplitude process oe slowly varying with 
respect to the local randan process. The properties of the dynamic systems can 
strongly influence the accuracy of the approximation, although the approxima- 
tion is accurate in the limiting case of slow variation of the amplitude pro- 
cess relative to the local process. For linear systems with low damping, the 
quasi-steady approximation is not valid in the nonlimiting case, since the 
values of the fourth-order flatness factor of the response are much closer to 
those of a Gaussian process than those of an amplitude modulated process. 

The response of dynamic systems to the sum of a product process and an 
independent Gaussian mean value process is examined. For linear systems the 
response can be developed by considering the response to the independent pro- 
duct and mean value processes separately. For nonlinear systems the coupled 
response can be developed by using the Fokker-Planck equation. 


Langley Research Center 
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APPENDIX A 


RELATIONS FOR CHARACTERISTIC FUNCTIONS 

Several relations which are used in the development of the characteristic 
function equations and the corresponding moment equations are derived in this 
appendix. The characteristic function is the Fourier transformation of the 
probability density function 


OO 

C(C) = E[e^^^] = I p(z) dz = F^p(z)^ 


(Al) 


The moments of a random process are obtained from the characteristic function 



C=0 



(A2) 


The development of the characteristic function equation from the Fokker- 
Planck equation requires several relations associated with the Fourier trans- 
formation. The first relation is the Fourier transformation of the derivatives 
of the probability density function 


F 



= i-^ C(C) 


(A3) 


The second relation is the Fourier transformation of the product of the inde- 
pendent variable and its probability density function 


f{z^ P(z)} = i~^ C(^) 

d^^ 


(A4) 


Both equations follow from the definition of the characteristic function 
(eq. (AT)). Another relation is obtained by combining equations (A3) and (A4) 


F 



= < 


dC (C) 
dC 


(A5) 
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The development of the moment equation requires the following relations 
for the derivatives of the products of the independent variable and the charac- 
teristic function: 



— c(o] = - 1)—— c(c) 

5=0 


c=o 


(A6a) 


(A6b) 


Equations (A6a) and (A6b) can also be expressed in terms of the moments of the 
random process by using equation (A2) 


— [C C(0] = fi-i®'-! e[z^~''] 

d?*” C=0 


' — C(0] = I {I - l)i^-2 e[z^-2] 

?=0 


(A7a) 


(A7b) 
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RESPONSE OF LINEAR DYNAMIC SYSTEMS - GENERAL R PROCESS 

The development in the body of the paper, which considers first-order 
filtering of white noise for the two components of the product process, can be 
extended to allow more general forms of the component processes. The relations 
for the response of linear systems to a product process whose R component 
process is generated by a general linear filtering of Gaussian white noise are 
developed in this appendix. This general form includes both the Dryden and, 
by use of rational approximations, the von Karman forms which are used to model 
atmospheric turbulence in aeronautical applications (ref. 10). The analytical 
development considers the response of a general linear system, both in coupled 
form for the development of the covariance matrix equation and in uncoupled 
form for the development of the general moment equation. 


Differential Equations 

The local R component of the product process is generated by a general- 
order linear filtering of a Gaussian white noise process. The corresponding 
state equation of the vector R process is 


rq = + da C2(t) (Bl) 

where A is the response coefficient matrix and d is the excitation coeffi- 
cient matrix of the state equation. The amplitude component is generated by 
a first-order linear filtering of an independent Gaussian white noise process 


s = -agS + ^2ag b ^-) (t) (B2) 

The excitation function of the dynamic system is the product of the amplitude 
ccmponent process and a linear combination of the ra variables, thus allowing 
a general form for the product process. The differential equation for the 
dynamic system is written as a linear state equation either in coupled or in 
uncoupled form 

4 

yj = Fjkyk + 9js<t>ara (B3a) 

4 

Zj = -MjZj + YjS<J)ai'a (No summation on j) (B3b) 

where 4> is a coefficient matrix which allows a general form for the contri- 
bution of the vector R process to the product process. The two sets of 
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coordinates are related by a linear transformation using the eigenvectors of 
the state coefficient matrix of the coupled system. 


Covariance Matrix Equation 

The differential equations for the canponents of the product process and 
for the coupled linear system are combined to form the state equation for the 
composite dynamic system. The conditional joint probability density function 
of this system satisfies the associated Fokker-Planck equation 
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— (sp) 
3s 


+ b^ 


3^p 
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3^P 
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- Fjk — (ykP) 
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3p 

3yj 


(B4) 


The corresponding characteristic function is the joint Fourier transform of the 
probability density function 


c(a,Pct,ilj |so,roio,yjo) = f{p(s, ra,yj | Sq, rao» yjo)) 


(B5) 


The equation for the characteristic function is the joint Fourier transform of 
the Fokker-Planck equation 
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3C 1 
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i9j^ j't'a 


32c 
3o 3 


(B6) 


The equation for the covariance matrix of the system response is obtained 
by taking the derivative of the characteristic function equation with respect 
to two of the transformation variables corresponding to the system response 
variables. The resulting covariance equation is 

3 

- — ^jS, + Fjk^kJl + F^k^kj = -9j<l)a%Jl - gjlttiaHaj (B7) 


where 


^kil = E[ykyjTj 

= E[srotys,] 
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The response terms on the left-hand side are identical to those in the corre- 
sponding equation for the response to Gaussian white noise excitation. The 
terms on the right-hand side are different, involving the H matrix of triple 
moments between the S process and one element of each of the R and 
Y arrays. The equation for these triple moments is obtained from the charac- 
teristic function equation by taking the derivatives with respect to the three 
transformation variables corresponding to the variables of the triple moments. 
The resulting equation is 

— - FjikHak = b^gj^tt^eRga (B8) 


where 

Rgot = E[rgrct] 


The covariance matrix for the vector R process is determined from the covari- 
ance matrix equation for the linear system of equation (Bl) 

3 

- ~ Rotg + AotYRyg + AgyRYOl ~ 
ot 

For the special case of a scalar R process, the preceding relations reduce 
to those developed in the body of the paper. The covariance matrix of the sys- 
tem response is obtained by solving the preceding relations in reverse order; 
equation (B9) gives the covariance matrix of the vector R process; equa- 
tion (B8) gives the H matrix of the triple moments; and equation (B7) gives 
the covariance matrix of the system response. 


Moment Equation 

The general moment equation for the composite dynamic system is developed 
by writing the equations for both the dynamic system and the vector R process 
in uncoupled form. The resulting composite dynamic system consists of the 
first-order differential equation for the amplitude component process 
(eq. (B2) ) , the set of equations for the uncoupled variables (q) of the 
vector R process, and the set of equations for the uncoupled linear system 

^ 2 ^^) (No summation on a) (BlOa)^ 

Zj = -PjZj + YjS^Jgqg (No summation on j) (BlOb) 


^The transverse velocity component of the Dryden spectral function 
requires special consideration since it introduces multiple eigenvalues. 
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where the -t, 6, and ^ arrays are obtained by transformation of the A, d, 

and 4) arrays of equations (Bl) and (B3b) . 


The joint probability density function of this composite dynamic system 
satisfies the associated Fokker-Planck equation 
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(BID 


The moment equation of the composite dynamic system is developed either from 
the Fokker-Planck equation or fran the corresponding equation for the charac- 
teristic function. The following notation is used for the joint statistical 
moments: 


Ejsl^ qi^l . . . 

The resulting moment 


z*| ^ « 2|.^ — M (k ; X"] 

equation for the composite dynamic system is 


(B12) 


X, n 


a=l j=1 


M (k;X-| , . . . ,Xg,;Vi , . . . ,Vn) 


— k(k — 1) agb2 m (k— 2;X'^,,.,,Xjj,; Vi ^ • r ^n^ 
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(B13) 
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The first term on the left-hand side of the moment equation corresponds to the 
transient response of the composite dynamic system. The remaining three sets 
of terms on the left-hand side correspond to the linear response terms of the 
conposite dynamic system. The first term on the right-hand side corresponds 
to the Gaussian white noise excitation term of the amplitude component process. 
The next two sets of terms on the right-hand side correspond to the white noise 
excitation terms of the vector Q process. These terms on the right-hand side 
pass the moments from lower to higher in the components of the excitation pro- 
cesses without changing the order in the system response variables. The final 
set of terms corresponds to the excitation term of the linear dynamic system. 
This set of terms passes the moments of given order in the response variables 
(and higher order in the amplitude and one component of the vector Q process) 
to the moments of higher order in the response variables. The moment equation 
can be solved numerically by the procedure described in the body of the paper 
for the case of a scalar R process. The resulting moments in the uncoupled 
response variables are transformed to those of the original dynamic system by 
using the linear transformation of the eigenvectors of the coefficient matrix 
of the original dynamic system. 
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Figure 3.- Fourth-order flatness factor of response of first-order system 

to product process. 




Figure 4.- Fourth-order flatness factor of first derivative of response of 

first-order system to product process. 



Figure 5.- Normalized fourth moment of response of second-order system 

to product process. 
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